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A  MATHEMATICAL  MODEL  OF  THE  ON-DECK 
HELICOPTER/SHIP  DYNAMIC  INTERFACE 


by 

J.  BLACKWELL  and  R.A.  FEIK 


SUMMARY 

A  mathematical  model  of  the  on-dock  helicopter/ship  dynamic  interface  has  been  developed  and  implemented  on 
an  ELXSI 6400  computer.  The  purpose  of  this  work  is  to  provide  a  capability  for  investigating  heUcopter/ship 
dynamic  interactions,  such  as  deck  clearances  on  landing,  swaying,  toppling  and  sliding  criteria,  and  tie-down 
loads.  Diiferont  helicopter  types  can  be  readily  examined,  given  their  undercarriage  representation.  The  model 
takes  account  of  arbitrary  ship  motion,  and  includes  features  such  as  brakes  on  orolf,  tyre  deformation  and  tyre 
sliding. 
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FIGURES 


NOTATION 


a  Acceleration  of  ship 

Cx,  Cy  Ty.  e  spring  coefficients  in  drag,  slide  directions 
C,,C2,  C3  Constants  used  in  Cx  and  Cy  (dependent  on  tyre  properties) 

Cz  Tyre  spring  constant  (zb  direction) 

d  Qleo  Deflection 

d0  Oleo  Break  Point 

d'  Tyre  compression  (zh  direction) 

F  Force  acting  on  helicopter  -  subscripts  tot  (totai),  g  (gravity),  u  (undercarriage), 
t  gnrust),  aer  (aerodynamic),  r  (restraining  devices) 

Fs,  Fb  etc.  Frames  of  reference 

F*  Undercarriage  force  without  sliding  or  dragging  included 

g  Acceleration  due  to  gravity 

g  Gravity  vector 

G  Moment  about  helicopter  eg.  Subscripts  as  for  F 
Gt,  G2  Damping  constants  for  oleo 

G*,  Gy,  Gz  Damping  constants  for  tyre 

h  Angular  momentum  of  helicopter  about  its  eg 

I  Matrix  determined  from  moments  of  Inertia  for  helicopter 

Ix,  Iy,  \  Moments  of  inertia  for  helicopter 

Iyy,  I„,  Iyj  Moment  of  inertia  crossed  terms  for  helicopter 

J'01  Matrix  required  for  conversion  from  <o' b  to  6'01 

K,  ,  K2  Spring  constants  for  oieo 

L  Oleo  load,  magnitude,  L, 

Lfc  Change  of  reference  frame  matrix  from  frame  s  to  b  (Appendix  A) 

Lx,  Ly,  Lj  Gear  moment  arms  (location  of  top  of  oleos  relative  to  helicopter  eg,  in  x,  y 

and  z  directions) 

L, ,  L2  Magnitude  of  oleo  load,  tyre  load 

M  Mass  ol  helicopter 

r  Position  vector  from  ship  eg  to  helicopter  eg 

r‘  Position  vector  from  helicopter  eg  to  tyre  contact  point  with  deck 
R  Reaction  torce  of  tyre  (normal  to  deck) 

Axa,  Aya  Tyre  deformation  in  fore/aft,  side  directions 

ii 


Aya,  Ay, 


Rate  of  tyre  deformation  in  fore/alt,  side  directions 


$>b,  0b,  \yb  Helicopter  attitude  (Euler  angles  from  earth  frame  to  helicopter  frame);  roll, 
pitch,  yaw 

<t>'01,  O'01,  v/cl  Relative  altitude  (Euler  angles  from  ship  frame  to  helicopter  frame);  roll,  pitch, 
yaw 

4>s,  0s,  Vs  Ship  ai  ,tude  (Euler  angles  from  earth  frame  to  ship  frame);  roll,  pitch,  yaw 

p*.  Uy  Drag  a  d  slide  tyre  friction  coefficients 

o >b  Aosolute  angular  velocity  ol  helicopter  (relative  to  earth) 

w'01  Angular  velocity  of  helicopter  relative  to  ship 

(os  Absolute  angular  velocity  of  ship  (relative  to  earth) 

SUBSCRIPTS  (AXES  SYSTEM) 

a  Aligned  axes  system 

b  Body  (helicopter)  axes  system 

e  earth  axes  system  (assumed  inertial) 

s  ship  axes  system 

SUPERSCRIPTS 

j  refers  to  j*1  undercarriage  gear,  eg.  Fu>  is  undercarriage  force  due  to  j1”  gear 
GENERAL 

Vectors  written  in  bold  type,  eg.  r 
Matrices  distinguished  by  a  tilde,  eg.  L 


1.  INTRODUCTION 


The  current  use  ol  helicopters  with  non  aviation  ships  creates  many  difficulties,  especially 
during  landing  in  the  presence  of  high  winds  and  rough  seas  (Ref.  1).  Some  of  these 
problems  are  being  partially  addressed  by  the  development  of  devices  such  as  the  "Lending 
Period  Designator"  (Ref.  2),  which  attempts  to  predict  relatively  quiescc.it  deck  motion 
intervals.  However,  once  the  helicopter  is  on  the  deck,  new  problems  arise.  The  deck  motion 
now  acts  directly  on  the  helicopter  landing  gear  and  may  cause  the  helicopter  to  osciilate  or 
slide,  making  deck  operations  hazardous.  In  high  sea  states,  the  helicopter  may  even  topple, 
causing  possible  damage  and  putting  deck  personnel  at  risk. 

The  Aircraft  Behaviour  Studies  -  Rotary  Wing  (ABS-RW)  Group  at  Aeronautical  Research 
Laboratory  (ARL)  is  currently  engaged  in  a  project  to  examine  the  on-deck  aspect  of  the 
helicopter/ship  dynamic  interface  problem.  This  document  describes  a  mathematical  model 
which  simulates  helicopter  interactions  with  a  moving  deck.  Rotor  thrust  and  other 
aerodynamic  forces  are  not  included  in  the  model  at  this  stage.  However  their  inclusion  is 
relatively  straightforward,  and  will  be  implemented  during  the  validation  stage. 

Section  2  briefly  outlines  the  problem  and  discusses  the  various  frames  of  reference 
utilised.  Sections  3  and  4  examine  the  forces  and  moments  which  act  on  the  helicopter  while 
on  an  arbitrarily  moving  deck.  These  comprise  mainly  the  undercarriage  forces,  which  are 
modelled  using  non-linear  equations.  Features  such  as  brakes  on  or  off,  tyre  deformation  and 
tyre  sliding  are  included  in  the  model.  In  Section  5,  the  linear  and  angular  accelerc  lions  are 
deduced.  Section  6  uses  the  mathematical  model  to  predict  helicopter  motions  under  a 
variety  of  deck  conditions.  Section  7  concludes  the  document,  whilst  Appendix  A  outlines 
the  unsteady  equations  req  lired  by  the  model. 

The  model,  which  is  not  specific  to  a  particular  helicopter,  was  implemented  on  an  ELXSI 
6400  computer.  A  standard  FORTRAN  version,  as  weif  as  a  simulation  language  version 
( ACSL ,  Ref.  3)  was  written. 


2.  OUTLINE  OF  PROBLEM 

Consider  a  helicopter  resting  on  a  ship.  Fomes  acting  on  the  helicopter  (eg.  gravity, 
undercarriage)  can  be  analysed  in  terms  of  forces  on  the  helicopter  eg  and  moments  about 
the  eg.  These  forces  and  moments  create  both  linear  and  angular  accelerations  of  the 
helicopter  relative  to  the  ship.  The  accelerations  can  be  integrated  numerically  to  determine 
position  and  attitude  of  the  helicopter  relative  to  the  ship  as  a  function  of  time,  for  various  ship 
motions.  Note  that  we  are  calculating  the  position  and  attitude  of  a  moving  reference  frame 
(helicopter)  relative  to  a  different  moving  reference  frame  (ship),  see  Fig.  1,  and  hence  the  full 
unsteady  equations  of  motion  are  required  (see  eg.  Ref.  4).  Those  unsteady  equations 
required  by  the  model  are  summarised  in  Appendix  A. 
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The  axes  systems  used  here  are  earth,  body  (helicopter)  and  ship.  They  are  rep  n.  rented  by 
subscripts  e,  b  and  s  respectively  (Fig  2).  A  further  axes  system,  represented  by  subscript  a, 
is  introduced  in  Section  3. 


3.  FORCES  ACTING  ON  HELICOPTER 

The  forces  acting  on  the  helicopter  are  : 

•  Weight  (gravity) 

•  Undercarriage 

•  Rotor  Thrust 

•  Aerodynamic  Forces 

•  Restraining  Devices 

The  riel  force,  Flot  is  the  sum  of  the  above  forces  : 

Ftoi  +  Fu  +  i  +  Fac  +  Ff  [11 

We  consider  here  only  the  we  jht  and  undercarriage  terms. 
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3.2  Undercarriage 

For  each  gear,  consider  the  model  shown  in  Fig.  3. 
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3.2.2  Side  and  Fore/Aft  Tyre  Forces 

All  the  undercarriage  wheels  are  assumed  to  point  in  the  direction  of  the  body  x-axis  at  this 
stage,  le.  no  provision  has  been  made  for  a  castoring  tail  wheel.  The  brakes  on  each  wheel 
can  be  either  full  on  (wheel  locked)  or  off  (wheel  can  rotate  freely). 


i)  Brakes  Full  On 

If  a  tyre  is  in  contact  with  the  deck  (L,  >  0),  any  movement  of  the  oleo  base  (wheel  axle) 
parallel  to  the  plane  of  the  deck  will  be  opposed  by  a  friction  force  at  the  point  of  contact 
between  the  tyre  and  the  deck,  which  will  deform  the  lyre  and  create  a  restoring  force.  When 
the  maximum  static  friction  force  is  exceeded,  the  wheel  will  skid  over  the  deck. 

It  is  necessary  to  define  a  new  reference  frame,  the  a-frame  (a  for  aligned),  obtained  by 
simply  rotating  the  ship-axes  frame  about  the  zs  axis  through  angle  t/el.  This  takes  account  of 

the  ship  and  helicopter  not  being  aligned  and  is  required  because  whilst  any  skidding  of  the 
helicopter  clearly  takes  place  in  the  plane  of  the  deck  (x,-ys  plane),  restoring  forces  are  likely  to 
differ  for  sliding  (skidding  in  the  direction  of  the  component  of  yb  in  the  xs-y5  plane,  the  side 
direction,  ya),  and  dragging  (skidding  in  the  direction  of  the  component  of  xb  in  the  xs-ys  plane, 
the  fore/aft  direction,  xa ).  See  Fig.  5.  The  representation  of  these  restoring  forces  for  the 
side  and  fore/alt  directions  is  a  function  of  tyre  geometry. 


The  L-matrix  from  frame  s  to  frame  a  is : 


/cosv/®1  siny™1  0\ 

-sm\/<'1  cosy'01  0 


(9) 


Given  the  updated  position  and  attitude  of  the  helicopter  during  a  time  history  calculation,  it 
is  possible  to  determine  the  position  of  each  oleo.  Assuming  that  the  tyre  base  is  at  its  original 
(ie.  previous  iteration)  position,  the  possible  tyre  deformation  in  the  plane  of  the  deck, 
(Axs,Ays)  is  deduced  and  hence  (Axa,Aya)  =  (cos\/0,Axs  +  sim/^Ay, ,  cos>/olAy5  -  siny'';lAxs)  is 

obtained  (Fig.  6).  If  however  the  tyre  should  leave  the  deck,  the  deformation  is  set  equal  to 
zero. 


Restoring  forces  in  the  xa  and  ya  directions  are  assumed  to  be  directly  proportional  to  tyre 
deformations  Axa  and  Aya  respectively  (Ref.  6),  and  are  given  by  (-CxAxa,  -CyAya).  The 
spring  coefficients  Cx  and  Cy  are  not  taken  as  constant,  but  are  assumed  to  vary  with  vertical 
tyre  deflection,  d‘  according  to  the  approximate  relation  (Loc.  Cit.) : 

Cx  -  C,  Md'  [10] 

C,  »  C2  -  C3d  [11] 


for  constants  C,,  C2  and  C3  which  depend  on  the  tyre  properties. 


The  damping  in  the  xa  and  ya  directions  is  assumed  to  be  directly  proportional  to  the  rate  of 
tyre  deformation  in  the  respective  directions,  Axa  and  Aya 
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The  total  undercarriage  forces  in  the  side  and  fore/alt  directions  may  also  contain  a 
component  of  the  oleo  load,  L,*  if  the  oleo  is  not  perpendicular  to  the  deck.  From  equations 
{8]  and  (9),  the  j5”  oleo  load  in  the  a-frame,  La'  is  given  by  : 

f  -cos<j>'olsinOrcl  \ 

L-i  =  Las  U1  =  simp'®1  L,<  [12] 

k  -cosb'clcosO'cl 

Thus,  the  total  force  acting  on  the  jlh  gear  (in  a-axes),  assuming  deck  contact  and  assuming 
no  sliding  or  dragging,  is  : 


-Cx*  Axa*  -  Gx*  Axa  -  cosb'^sinO'01  L,1 
"CyAyJ  -  Gy*  Aya  +  sinb^'L,1 
-cos<i>'olcosQ'el  L,’ 


total  force  in , ore/aft  direction 
total  force  in  s  de  direction  [13] 
'otal  lorce  normal  to  deck 


if  the  j*1  tyre  is  off  the  deck,  Fuj  is  clearly  zero.  The  foices  for  a  sliding  or  jragging  tyre  are 
discussed  next. 

The  limiting  friction  forces  in  the  fore/aft  and  side  directions  are  p,1  i  ’J  jnd  p  j  R'  respectively. 
R',  the  force  normal  to  the  deck,  is  given  by  the  modulus  of  'Equation  13),  and  is  equal 
to  cos4»ro1  cosG'"1  L,i .  If  the  modulus  of  the  fore/aft  force  (or  side  force)  n  t-qua  non  [13]  exceeds 
the  limiting  drag  (slide)  friction,  then  the  fore/aft  (side)  force  is  set  equ.T  to  fl  >r.  I  mi  ting  friction 
and  the  tyre  base  is  moved  rigidly  with  the  oleo,  such  that  the  d  ;<  >rrnatir  n  (ii  (he  plane  of  the 
deck)  of  the  previous  iteration  is  retained.  Thus,  no  new  tyie  dt.crmalion  (iV  (he  plane  of  the 
deck)  is  assumed  to  occur  and  the  tyre  drags  (slides).  It  ’he  f pre/aft  (side)  force  does  not 
exceed  the  limiting  value,  the  tyre  base  is  not  moved  and  instead  the  tyro  deforms.  See 
Fig.  7. 

Final  equations  for  the  jth  gear  load,  F^1  are : 


where : 

RJ  =  cos^cosO™!,*  [15] 


|Fua'  |x  |F„af  jx  |  <  px‘  R1  (tyre  doesn't  drag) 

jnxi R1  .  sign(F*al  |x)  |F„a»  |x  |  >  px*  R>  (tyre  drags) 
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117] 


K'ly  Ka’U^Hy'R' 

03  ly  =  juy'R1  •  Sign(Fu0*  |y)  \K>  ly  |  >  Hy*  R' 

with : 

K*'lx  =  -cxf  Axa'  -  Gx'  Axa  -  cos<j.'0,sinGrel  L,» 
F^'ly  =  "Cy  Aya*  •  Gy*  Aya  +  sin4>,el  L,* 
sig»(F)  =  j^j- 

unless  the  ith  gear  is  off  the  ground,  in  which  case  F„  *  =  0. 


(tyre  doesn't  slide) 
(tyre  slides) 


[18] 

[19] 

[20] 


ii)  Brakes  Off 

For  the  case  of  brakes  off,  the  undercarriage  wheel  (assuming  frictionless  bearings  and 
neglecting  rolling  friction)  can  rotate  freely  in  the  xa  direction  and  so  does  not  drag  or  deform. 
This  affects  the  xa  forces  as  shown  below.  The  sliding  or  deformation  in  the  ya  direction  is 
assumed  to  oe  unaffected  by  the  status  of  the  brakes,  and  so  behaves  as  in  the  preceding 
section. 

Given  ttv.  updated  position  of  the  oleo,  the  tyre  is  assumed  to  move  in  the  xa  direction  such 

that  it  is  it  iine  with  the  oleo,  ie.  no  deformation  takes  place,  and  hence  no  restoring  force  acts 
in  the  xa  erection.  Thus,  the  xa  component  of  FUa*  contains  only  the  xa  component  of  oleo 

load  anc.  ;3  modified  from  equation  [16]  to  read  : 

F^'l,  =  -  cos^sinO'^L,)  [21] 

with  y,  ana  za  components  assumed  to  be  unaffected. 


Fer  brakes  on  or  off,  the  oleo  load  L,*  is  determined  from  the  oleo  and  tyre  zb  deflections. 

Those  j.e  obtained  from  the  helicopter  position  and  attitude,  as  are  the  tyre  displacements  in 
the  plane  of  the  deck,  Ax.,1  and  AyaJ.  Summation  of  Fj  over  the  separate  gears  leads  to  the 
resultant  undercarriage  lorce,  Fu : 


F  J  -  I 
rus  -  Lsa 


=  X'v 


.  Pva> 

[22] 

[23] 

4.  MOMENTS  ACTING  ON  HELICOPTER 

The  forces  discussed  in  Section  3  do  not  in  general  act  at  the  helicopter  eg,  and  hence 
moments  are  created  about  the  eg.  The  net  moment  about  the  eg,  G10t  is  defined  as  the  sum 

of  the  moments  due  to  the  afore  mentioned  forces  : 
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Gg  +  Gu  +  °i  +  G,«  +  Gr 


(24) 


G 


toi 


with  only  the  weight  and  undercarriage  terms  considered  here. 


4.1  Weight 

The  moment  due  to  gravity,  Gg,  is  zero  since  gravity  acts  at  the  helicopter  eg. 


4.2  Undercarriage 

The  net  force  acting  on  the  jth  gear  is  given  by  F J ,  equations  [14]  -  [17].  The  moment 
about  the  helicopter  eg  due  to  this  force,  Gj  is  : 

G0j‘  =  x  r'f  [25] 

where  r'i  is  the  position  vector  from  the  helicopter  eg  to  the  jUl  tyre  contact  point  with  the  deck 
(Fig.  8). 

Summation  of  Gu*  over  the  separate  gears  leads  to  the  resultant  moment  about  the 
helicopter  eg  due  to  undercarriage  forces,  Gu : 


GUb'  ---  L*  GuJ 

[26] 

G%  =  XGub 

[27] 

i 


Given  the  force  acting  on  the  helicopter  eg  and  the  moment  tending  to  rotate  the  helicopter 
about  its  eg,  the  linear  and  angular  accelerations  can  be  determined.  Numerical  integration 
(fourth-order  Runge-Kuita)  is  then  used  to  deduce  the  position,  attitude  and  velocity  of  the 
helicopter  with  respect  to  the  ship. 


5.  EXPRESSIONS  FOR  LINEAR  AND  ANGULAR  ACCELERATION 

Not  only  is  the  position  and  attitude  of  the  helicopter  relative  to  the  ship  changing  with  time, 
but  the  ship's  position  and  attitude  relative  to  the  earth  is  also  constantly  changing.  Hence  it  is 
necessary  to  use  the  full  unsteady  equations  of  motion  which  take  account  of  the  moving 
reference  frames. 


5.1  Linear  Acceleration 

If  r  is  the  position  vector  of  the  helicopter  eg  relative  to  the  ship  eg  (Fig.  9),  then  from  Ref.  4  : 

“  I* lo!«  •  * 

r>  “  ~  -  -  2»* x r,  -  x (»* x g  [28] 
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where : 


a  =  acceleration  ot  ship  eg  (relative  to  earth) 

«s  =  angular  velocity  of  ship  (relative  to  earth) 

M  =  mass  of  helicopter 
Ftot  =  total  force  acting  on  helicopter 

Thus,  given  the  net  force  acting  on  the  helicopter  eg  in  ship  axes  (calculated  in  Section  3),  the 
ship's  acceleration,  angular  velocity  and  angular  acceleration  relative  to  earth  (input  to  model), 
and  the  current  position  and  velocity  of  the  helicopter  eg  relative  to  tne  ship  (deduced  by 
numerical  integration  over  the  previous  time  step),  we  can  calculate  the  acceleration  of  the 
helicopter  relative  to  the  ship,  r, .  Numericai  integration  is  then  used  to  determine  updated 
values  for  the  velocity  and  position  of  the  helicopter,  rs  and  rs . 


5.2  Angular  Acceleration 

If  wb  is  the  absolute  angular  velocity  of  the  helicopter  (relative  to  earth),  and  Gt0,  is  the  net 
moment  about  the  helicopter  eg,  calculated  in  Section  4,  then  the  angular  acceleration  of  the 
helicopter,  <ob  is  given  by  (see  Appendix  A) : 

IGrw,,  -  o£x(Icob)]  [29) 


where  I  is  obtained  from  the  helicopter  moments  of  inertia,  and  is  given  by  : 


I  = 


(  I*  ’^xy  ”Vz 

*Vy  V  'Vi 

V,  "Vz  'Vz  V  ) 


[30] 


Numerical  integration  is  then  used  to  determine  an  updated  value  for  o>b,  which  is  used  to 
deduce  the  angular  velocity  of  the  helicopter  relative  to  the  ship,  a/01  from  the  relation  : 


=  %  '  Lbs  “1 


[31] 


The  Euler  angle  rates  are  developed  in  Ref.  4,  and  as  shown  in  Appendix  A  are  given  by  • 


fSl 


Kr 


J'01  co'g1 


[32] 


numerical  integration  can  now  be  used  to  determine  <{>'cl,  0'el  and  y™1,  the  attitude  of  the 
helicopter  relative  to  the  ship. 


Given  updated  values  of  relative  position  and  attitude,  the  whole  procedure  is  now 
repeated. 
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6.  APPLICATIONS 


Consider  a  helicopter  "dropped”  from  a  horizontal  position,  onto  a  ship  deck  which  has  a 
variety  of  attitudes  and  motions.  For  the  purpose  of  illustration,  we  restrict  ourselves  to 
</“'  =  0°  (helicopter  and  ship  aligned),  and  the  ship  is  allowed  to  deviate  from  the  horizontal 
position  in  roll  only  (Fig.  10).  The  helicopter  considered  has  two  main  landing  gears  with 
brakes  ON  ,  and  a  single  tail  gear  without  brakes  (ie.  tail  brake  taken  as  OFF  in  the 
mathematical  model).  Rotor  thrust  and  other  aerodynamic  forces  are  taken  as  zero  'or  this 
idealised  case. 

The  following  cases  aie  considered  : 

case  1  -  deck  horizontal  and  stationary 

case  2  -  deck  at  10°  roll  and  stationary 

case3  -  deck  rolling  sinusoidally  from  -10°  to  +10°  with  period  10  seconds 
The  various  parameters  representing  the  helicopter  are  listed  in  Tables  1  and  2. 


Table  1.  Helicopter  Parameters 


Parameter 

SI 

Imperial 

M 

9.0  xIO3  Kg 

2.0  xio4  lb 

I* 

9.0  x  1 03  Kg  mz 

2.1  xio5  lb  ft2 

Jy 

6.3  xIO4  Kgm2 

1.5  x10s  lb  ft2 

k 

6.0  xIO4  Kgm2 

1.4x10®  lb  ft2 

2.5  xIO3  Kgm2 

5.9x10"  lb  ft2 

^xy=^yx 

0.0  Kgm2 

0.0  lb  ft2 

0.0  Kgm2 

0.0  lb  ft2 

The  mathematical  model  developed  earlier  is  used  to  predict  the  variation  of  oleo  load  and 
oleo  deflection  with  time  for  the  above  three  cases. 

Case  1  -  See  Fig.  1 1 .  The  oleo  load  and  deflection  variation  for  gears  1  and  2  are  the  same 
due  to  symmetry.  The  helicopter  comes  to  rest  in  about  three  seconds,  with  greater 
compression  of  the  tail  wheel  oleo  than  the  main  gear  oleos. 

Case  2  -  See  Fig.  12.  Due  to  the  initial  horizontal  attitude  of  the  helicopter,  the  left  main 
gear  (1)  makes  contact  with  the  deck  first,  followed  by  the  tail  gear  (3)  and  finally  the  right  main 
gear  (2).  After  making  deck  contact,  gear  1  actually  leaves  the  deck  for  a  short  while.  Case  2 
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takes  longer  to  damp  out  than  case  1 ,  as  expected,  and  the  helicopter  comes  to  rest  after 
about  six  seconds,  with  the  right  main  gear  (2)  more  compressed  than  the  left  main  gear  (1). 

Case  3  -  See  Fig.  13.  The  relatively  gentle  ship  rolling  motion  has  little  effect  on  oleo  load 
and  deflection  during  the  initial  impact  (compare  Figs.  12  and  13),  however  by  five  seconds, 
when  the  ship  deck  is  now  inclined  the  other  way,  a  reversal  of  the  loads  on  the  main  gears 
has  taken  place.  If  the  variation  of  oleo  load  and  deflection  is  examined  over  a  greater  time 
span  (Fig.  14),  it  is  seen  that  a  periodicity  in  the  helicopter  motion  soon  develops,  with  period 
10  seconds  equal  to  that  of  the  ship.  The  helicopter  Is  seen  to  oscillate  between  the  two  main 
gears,  with  the  tail  gear  oleo  load  and  deflection  remaining  virtually  constant. 


Table  2.  Helicopter  Landing  Gear  Parameters 


Parametei 

MAIN  GEAR 

TAIL  GEAR 

K» 

3.6  x  10s 

Nrrv2 

(7.5  x103  lbf/ft2) 

1.8  x  10s  Nm'2 

(3.8  xiO3 

lbf/ft2) 

k2 

5.0  x  10® 

Nnv2 

(1.0  x10s  lbf/ft2) 

2.0  x  10®  Nm  2 

O 

»— 

X 

CM 

lbf/ft2) 

Gi 

O 

X 

CM 

CM 

Nsnr' 

(1.5  xIO3  Ibf  s/ft) 

1.8  xio4  Nsnr1 

(1.2x103 

Ibf  s/ft) 

°2 

O 

X 

o 

Nsnr1 

(2.7  X103  Ibf  s/ft) 

3.0  XIO4  Nsm'1 

(2.1  x  103 

Ibf  s/ft) 

do 

0.238 

m 

(0.781  ft) 

0.335  m 

(1.099 

ft) 

c2 

6.0  x10s 

Nor1 

(4.1  xIO4  lbf/ft) 

8.0  x10s  Nm-' 

(5.5  xIO4 

lbf/ft) 

c, 

1.4  x  10® 

Nm"’ 33 

(6.5  xIO4  lbf/ft'33) 

8.0  x10s  Nm-’-33 

(3.9  xIO4 

lbf/ft’33) 

^2 

4.5  x  10s 

Nm"' 

(3.1  xIO4  lbf/ft) 

6.0  x10s  Nm-' 

(4.1  x  104 

lbf/ft) 

c3 

1.2  x  10® 

Nm2 

(2.5  x  1 04  lbf/ft2) 

1.0x10®  Nnr2 

(2.1  xIO4 

lbf/ft2) 

2.0  x  10" 

Nsnr1 

(1.4  xIO3  Ibf  s/ft) 

2.0  xio4  Nsm-’ 

(1.4x103 

Ibf  s/lt) 

GV 

2.0  x  104 

Nsnr’ 

(1.4x103  Ibf  s/ft) 

2.0  xIO4  Nsnr’ 

(1.4  XIO3 

Ibf  s/ft) 

G* 

1.0  x  104 

Nsm-' 

(685  Ibf  s/ft) 

1.0  xIO4  Nsm-’ 

(685 

Ibf  s/ft) 

Lx 

1.5 

m 

(4.92  ft) 

3.3  m 

(10.83 

«) 

Ly 

1.4 

m 

(4.59  ft) 

0.0  m 

(0.0 

ft) 

l2 

0.8 

m 

(2.62  ft) 

0.6  m 

(1.97 

M) 

P* 

0.7 

(0.7) 

0.7 

(0.7) 

0.7 

(0.7) 

0.7 

(0.7) 

A  possible  further  application  of  the  model  is  to  examine  the  varying  distance  between  the 
ship  deck  and  a  part  of  the  helicopter.  For  instance,  the  Sikorsky  SH-60B  Sea  Hawk  has  a 
radome  located  under  its  chin  which  could  possibly  come  close  to  the  deck  during  a  heavy 
landing,  or  while  the  helicopter  is  on  the  deck  during  rough  seas.  In  Fig.  15,  the  distance 
between  the  deck  and  a  hypothetical  radome  (location  [3.5,  0, 1 .5jm  relative  to  the  helicopter 
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eg),  is  plotted  for  the  conditions  of  case  3.  More  realistic  situations  would  require  the  use  of 
measured  ship  motions  as  forcing  inputs. 


7.  CONCLUDING  REMARKS 

A  mathematical  model  has  been  described  which  models  a  helicopter  on  an  arbitrarily 
moving  deck.  The  model  has  been  successfully  implemented  on  an  ELXSI 6400  Computer. 
Results  so  far,  for  a  deck  at  rest  and  oscillating  sinusoidally,  agree  with  expected  trends.  It  is 
planned  to  valioate  the  model  using  actual  helicopter  data  in  the  near  future. 

The  oleo  load  and  deflection  have  been  examined  here,  but  other  variables  of  interest, 
such  as  deck  clearance  of  a  portion  of  the  helicopter,  or  clearances  with  nearby  obstacles  di  >•' 
to  swaying,  can  be  easily  determined.  The  model  can  also  b6  used  for  developing  sliding  and 
toppling  criteria. 

The  addition  of  restraining  devices  (ti**  jown  and  RAST*  loads)  is  desired  to  extend  the 
capabilities  of  the  on-deck  mode;,  it  is  also  intended  to  include  rotor  thrust  and  aerodynamic 
forces  in  the  model,  and  a  stability  augmentation  system  (SAS),  so  that  actual  landings  on  an 
arbitrarily  moving  deck  can  be  studied.  Other  assumptions,  for  example  that  oteo  loads  act  in 
the  direction  of  the  body  2-axis,  can  be  relatively  easily  relaxed. 

The  model  is  readily  adaptable  to  different  helicopter  types.  In  particular,  the  undercarnage 
representation  can  be  varied,  either  by  changing  the  parameters  in  the  present  model,  or  by 
providing  an  alternative  representation,  without  changes  being  required  in  other  parts  of  the 
model. 

Finally,  the  model  has  been  developed  so  as  to  facilitate  the  incorporation  of  either  part  or  all 
of  it  into  wider  system  studies  of  the  helicopter/ship  dynamic  interface. 
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APPENDIX  A 


Unsteady  Equations  of  Motion 

Due  to  the  non-steady  nature  ol  the  on-deck  problem,  the  unsteady  equations  of  motion 
are  required.  For  a  full  discussion  of  these  equations,  see  eg.  Ref.4.  The  main  relationships 
used  in  this  document  are  listed  below,  with  the  relevant  section  number  from  Ref.  4  listed  in 
brackets  where  appropriate.  The  earth  frame  is  assumed  to  be  inertial  throughout. 

Position  and  Velocity  Vectors  in  Different  Reference  frames 

Consider  two  reference  frames  eg.  b  and  s.  If  rb  is  a  vector  viewed  from  frame  b,  and  rs  is  the 
same  vector  viewed  from  frame  s,  then  rb  and  r,  are  related  to  one  another  as  follows : 

ib  =  ^  rs  (§4.4)  [A1] 

where  LbS  is  the  change  of  reference  frame  matrix  from  s  to  b  given  by : 

/  cosO'£flcos'/°l  cosGrolsinyro1  -sin9,<H 

Lb,  =  sin^sinO^cos^-cos^siny01  sin</olsinG'olsiny,ol+cosy,olcosi)>r<!l  sin^'cosO'01 

^  cos4>,Qls  in  Orelcos\}/0,-i-si  cosd'0lsinG,0lsiny'0l-sinci>,Glcosi{;r01  cos0ro,cos<j>">l 

and  (j)'®1,  O'01  and  i/01  are  the  Euler  angles  that  take  frame  s  to  frame  b. 

Note  that : 

Lsb  =  LbTs  (§4.4)  [A2] 

The  rate  of  change  of  vector  r  when  viewed  from  frame  b,  rb  is  related  to  the  rate  of  change 
of  the  same  vector  when  viewed  from  frame  s,  r,  by : 

rb  =  Lt»  r,  -  <  x  rb  (§4.6)  [A3] 

where  c o'01  is  the  angular  velocity  of  frame  b  relative  to  frame  s. 

Angular  Rates  In  Different  Reference  frames 

If  G  is  the  moment  and  h  the  associated  angular  momentum  of  a  body  about  its  eg,  then  in 
body  axes,  b : 

Gb  =  hb  +  x  hb  (§5.6)  [A4] 

providing  that  the  origin  of  the  body  axes  is  the  body  eg.  <ob  is  the  absolute  angular  velocity  of 
frame  b  (relative  to  earth). 

For  a  rigid  body,  the  angular  momentum  is  given  by . 

h0  =  R  (§5.4)  [A5] 

where  I  is  obtained  from  the  body  (helicopter)  moments  of  inertia : 


f  I,  ~Vy  '^xz  s 

I  =  -Ixy  !»  -lyz  IA6] 

-*yi  *2  J 

and : 

Ix  =  J(y2  +  z2)dm  etc.  |A7) 

•xy  =  }xydm  etc.  [A8J 

Equations  [A4]  and  [A5]  result  in  : 

Gb  =  I  <i)b  +  c<£  *  (I  )  (§5.6)  [A9] 

Finally,  the  rate  of  change  of  the  relative  Euler  angles  are  related  to  the  relative  angular 

velocity  by : 

( <j,r°'  ^ 

6'°'  =  >'  <  (§5.2)  [A10] 

wv 

where  : 

/■  1  sinij>‘oltan0ro1  cos^tanO'01  \ 

j,el  =  0  cos<f>ro*  -sin$rol  [AH] 

x  0  sin<{>'olsccO'01  cos^'sccO'01  J 

and : 

<  =  (Ob  -  Ltj  (0*  (§5.2)  [A12] 

with  <os  being  the  absolute  angular  velocity  of  frame  s  (relative  to  earth). 
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L,  =  Oleo  Load 

Lj  =  Tyie  Load 

d  =  Oleo  Compression  (zb  Direction) 

d’  «  Tyre  Compression  (zb  Direction) 

d+d*  =  Total  Gear  Deflection  (Oleo  +  Tyre) 


Figure  4.  Oleo  and  Tyre  Loads  and  Deflections  :n  Body  z  Direction 


Line  of  Oleo 


Figure  6.  a)  Tyre  Deformation  in  Side  (ya)  Direction,  b)  Tyre  Deformation 
in  Fore-Aft  (x  )  Direction 


Fu^ly  >  HyJ  RJ 
Tyre  slides 

No  additional  deformation 


Filly 

No  sliding 

additional  tyre  deformation 


Figure  7.  a)  Tyre  Sliding  and  b)  Tyre  Deformation  in  Sideways  Direction  (Tyre 
Dragging  and  Deformation  in  Fore-Aft  Direction  is  Similar) 


Figure  10.  Helicopter  "dropped"  onto  a  Rolling  Ship  Deck 


Figure  11.  Oleo  Deflection  and  Load  for  Helicopter  "Droppers"  onto  Ship  Deck. 
Deck  Horizontal  and  Stationary 


Figure  13.  Oleo 
Deck  Rollir 


Tlmo 

c=> 


flection  and  Load  for  Helicopter  "Dropped"  onto  Ship  Deck, 
linusoidaily  from  -10°  to  +10°  with  Period  10  Seconds 


Tlmo 

(•) 


Figure  14.  ^j^ng  j  Ship^Deck. 


Figure  15.  Deck  Clearance  of  Radome. 

Deck  Rolling  Sinusoidally  from  -10°  to  +10°  with  Period  10  Seconds 
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